%%
circle1 = [
    
%%
imshow(i);
x = edge(i,'canny',[],2);
a = double(x);
figure;
imshow(x);
figure;
imshow(a);

smallR = 6;
bigR = 8;

small = fspecial('disk', smallR);
big = fspecial('disk', bigR);

smallPadding = ((bigR - smallR)*2) / 2;

for j=1:2
   small = [zeros(size(small,1),smallPadding)';small'];
end
for j=1:2
   small = [small';zeros(size(small,1),smallPadding)'];
end

mask = big - small;
imshow(big - small,[]);

Imout = convolve2(a,mask,'same'); 
figure;
imshow(Imout,[]);

[r,c,v] = findpeaks(Imout); % add thresh for r & c 

ok = v > 0.2; % play around to change thresh
r = r(ok);
c = c(ok);

hold on;
plot(c,r,'g*');


